Deconfinement energy threshold: analysis of hadron yields at 11.6 A GeV 
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We analyze the hadron yields obtained at the AGS in the range 11-11.6 A GeV and discuss 
strategies to identify possible deconfinement at this energy scale. These include consideration of 
chemical non-equilibrium at hadronization, and the study of (multi)strange hadrons. We find that 
the totality of experimental results available favors the interpretation as hadron freeze-out at the 
phase boundary between confined and deconfined phase. 
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One of the most interesting issues, in the field of rel- 
ativistic heavy ion collisions, is the understanding of 
the thresholds in reaction energy, and system size, be- 
yond which the formation of the color deconfined par- 
tonic state occurs dominantly. In the central collisions 
of Au-Au ions at BNL-RHIC, at the top energy range 
(v' s nn = 200 GeV), it is generally believed that a color 
deconfined state has been formed Q. The strange an- 
tibaryon production systematics lead us to believe that 
this is also the case at the top energy in central Pb-Pb 
ion reactions at y/s^ = 17.2 GeV (158 A GcV/c Pb 
beam colliding with fixed target) at CERN-SPS 0. Ex- 
perimental exploration of lower energy collisions at SPS 
down to ^/snn = 6.3 GeV (20 A GeV/c) suggests a pos- 
sible change in strangeness production pattern In 
this analysis, we address the experimental program at 
BNL-AGS at = 4.8 GeV (11.6 A GeV/c). Our 

objective is to apply the methods we used in the study of 
the SPS and RHIC data in order to identify similarities 
and differences in the hadron production pattern. 

The tool used in this study of hadron production is 
the Statistical Hadronization (SH) model introduced by 
Fermi in 1950 5]. In 50 years, SH has matured to a 
full fledged tool, in the study of soft strongly interacting 
particle production, capable to describe in detail hadron 
abundances once the statistical grand canonical method 
and the full spectrum of resonances is included [6]. The 
key SH parameters of interest are the temperature T and 
the baryochemical potential /ib- It is generally accepted 
that as the energy of the colliding nuclei varies, the prop- 
erties of the hadronic gas (HG) and quark-gluon plasma 
phases (QGP) are explored in a wide domain of T and 
Hb , see Fig. 2. in @. 

Two values of temperature can be determined in sta- 
tistical hadronization studies: chemical freeze-out T c h, 
which determines particle abundances, and a thermal 
freeze-out T t h, defined by the condition where the mo- 
mentum spectra stop evolving, i.e., particles stop inter- 
acting elastically. We analyze here solely the yields of 
hadronic and in particular strange particles and, hence 
from now, on we imply T — T c h. For T to be uniquely de- 
termined for all particles considered the chemical freeze- 



out must occur rather fast. This requires that particles 
stop interacting just after formation (hadronization). In 
fact, the same freeze-out condition explains, at high en- 
ergy SPS 0, and at RHIC @, both particle abundances 
and spectra. This is the so called single freeze-out model 
expected to apply in presence of sudden hadronization of 
a rapidly expanding supercooled quark-gluon fireball [t| . 

The single freeze-out model, allowing for chemical non- 
equilibrium, is consistent with the observation by invari- 
ant mass method of abundantly produced hadron res- 
onances by, e.g., the STAR and NA49 collaborations 
|ToL ITU . This results demonstrates that resonance decay 
products in essence have not scattered before their ob- 
servation [24j . Hence resonance decay spectra contribute 
to the stable particle spectra, altering the purely ther- 
mal shape. An analysis of this modified shape yields 
the universal freeze-out condition 0, On the other 
hand, studies in which resonance decays are for simplic- 
ity ignored have as result a different thermal freeze-out 
condition for practically each particle considered, which 
differs from the chemical freeze-out condition obtained 
using hadron yields. These may serve as a comparison 
benchmark but cannot be used otherwise. 

The question of chemical equilibration arises in the 
context of chemical freeze-out studies. The introduction, 
into the data analysis, of the strangeness phase space oc- 
cupancy parameter "f s was motivated by the recognition 
that strangeness (or cquivalcntly strange hadron) build- 
up in microscopic reactions will rarely lead to abundance 
expected in chemical equilibrium, j s = 1 12] . The pa- 
rameter 7 S characterize, independently of temperature, 
the particle-antiparticle pair yield. Subsequently, the 
ligh t quark phase space occupancy j q was introduced 
|13| . This allowed to describe particle yields for the case 
of fast hadronization. Namely, the phase space density of 
quarks in QGP is very different from that obtained eval- 
uating the yields of valance quark content in final state 
hadrons. 

Thus, even if during its temporal evolution a QGP fire- 
ball converges to chemical equilibrium, the final hadron 
state generated on the scale of a few fm/c should emerge 
showing a pattern of subtle deviations from chemical 
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equilibrium which contain interesting information about 
the physics we are exploring |14Ul5| . The claim that, in 
general, overall chemical equilibrium is reached in rela- 
tivistic heavy ion collision is based on an analysis that 
presupposes this result [l6|. In this paper, we show 
that while chemical equilibrium hypothesis for the 11 A 
GeV reactions can be considered, the experimental re- 
sults available today favor chemical under saturation of 
both light and strange quark phase spaces, with both 

7s, Iq < I- 

Let us first explain how the different parameters of 
the SH model can be determined using the experimental 
data in a step by step process, rather than in a global fit, 
which erases much of the physics insight. We will show 
that it is possible to describe a partial subset of data 
with a partial subset of parameters since certain types 
of particle ratios are sensitive to subsets of statistical 
parameters. Of particular importance in this discussion 
is that the chemical non-equilibrium parameters 7s,7q, 
require to be determined ratios of certain rarely produced 
particles. 

a) The chemical fugacities X q , A s , see \v^, an d A/ 
(see e.g Eq. (^below) can be determined rather precisely 
and independent of other statistical parameters consid- 
ering the anti-particle to particle ratios. There are many 
such ratios available allowing to fix the three parameters, 
check consistency and to make ratio predictions based on 
a partial data set. Success of this part of SH model anal- 
ysis has no predictive power regarding the value of the re- 
maining 4 statistical parameters. Alas, we see again and 
again the argument to the contrary, in particular some 
workers claim based on the success of describing particle 
to antiparticle ratios that there is chemical equilibrium, 
7s = 7g = 1, which is plainly wrong. 

b) If one considers from that point on only the prod- 
uct of particle with antiparticle yield, to a very good ap- 
proximation one obtains reduced yield ratios which are 
independent of A ? , A s and Aj . Than to determine the 
ratio ^s/lq (or equivalently 7„ when j q = 1) we consider 
ratios of particles with unequal strangeness content, e.g. 
0/VK+K - oc 7 s /7 9 /i(T) In such a ratio generally there 
is a correlation between 7s /7g and T and thus at least 
two such ratios are required to determine "f s /lq and T, 
see for example Ref . [l7j , figure 4 . 

c) Similarly comparing yields of particles with differ- 
ent quark number content, e.g. mesons with baryons 
one identifies the value of 7 g , as noted in e.g., Section 
IVD in Ref. 0. Consider here as example AA/K+K~ oc 
7<?/2(T); we see that at least a second such ratio is re- 
quired, e.g. \/s - S + /0 or/and NN /ir + n~ etc. to deter- 
mine both 7 g and T or/and we take T from the study of 
7s/7 9 and T above. We expressly did not introduce above 
pion yields as these derive from resonance decays includ- 
ing heavy baryons, and thus their yield is not strictly the 
yield of mesons. Moreover the pion yields are theoreti- 
cally most uncertain considering possible extension of the 
hadron mass spectrum to higher mass, and uncertainties 
about number of pions produced in cascading decays of 



heavy resonances. In absence of experimental data re- 
garding multistrange particles there is little choice but 
to compare the yields of pions to nucleons in the analysis 
of 7 9 and T. 

d) Each individual total particle yield is proportional 
to the volume parameter V aside of a strong dependence 
on T, and a lesser dependence on all the other SH pa- 
rameters. This allows to fix V and test the consistency of 
the individual findings about the other parameters made 
above. 

In a more modern approach one makes a global and 
simultaneous fit of all parameters minimizing the global 
error (x2). This can only work if there is sufficiently 
wide scope of experimental data as required in the indi- 
vidual steps above. Practical experience shows that an 
incomplete data set combined with the presence of a sig- 
nificant measurement error makes the task of unfolding 
7s, 7 9 and T difficult. In this environment some workers 
think that it is better to assume j q = 1, we will discuss 
this at length below. Here we note that the predictive 
power of SH model can be tested: we can create choos- 
ing a set of SH model parameters a set of particle yields, 
and add in random errors. If the requirements on avail- 
ability of measurements identified in a)-d) are respected 
in the data set considered, one can fit such an artificial 
data set to obtain the 'creating' set of statistical param- 
eters. However, when large yield errors are introduced, 
and key particles eliminated, the expected solution will 
not always be found in the likelihood analysis. 

Considering such studies it is important to remember 
that the SH model is not being verified or tested, its 
validity can be seen in its predictive power which encom- 
passes particle yields that vary by typically four orders 
of magnitude, varying between C(/.oo) and C(oo//). The 
issue before us is, if and when we can use statistical sig- 
nificance of the fit to infer the values of the parameters, 
and the insights obtained in such model studies are af- 
firming this. For this reason it makes sense to look for 
most likely solutions using a full parameter set motivated 
by theoretical considerations. 

The parameter "f q has been introduced into SH model 
not as a fit parameter, but in consequence of a theoret- 
ical development. It is the quantity which is expected 
to differ from unity considering hadron freeze-out from 
a rapidly expanding system and/or undergoing a sharp 
phase transformation. In such consideration 7 9 = 1 only 
when there is a relatively long time available for chemi- 
cal re-equilibration. Moreover, a large value of "f q allows 
the hadron gas phase to absorb without volume increase 
the enhanced entropy content of the deconfined phase in 
which color bonds are broken. The ability of the two 
models (j q = 1 and j q ^ 1) to describe data should be 
compared in a quantitative way in the study of statisti- 
cal significance. Evidence for j q > 1 constitutes evidence 
for a rapidly evolving system emerging from a deconfined 
phase PS |3 . 

Statistical significance |21|. defined as the probability 
of the data fit to a model to be of the obtained quality, 
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provided that the model under consideration is 'true', 
and errors are purely due to experiment, is an appro- 
priate tool for the comparison of model variants. This 
approach takes into account both the number of fit pa- 
rameters and the 'goodness' (%2) of the fit. In the ab- 
sence of the experimental data necessary for strict falsi- 
fication, it has been argued by others that a comparison 
of statistical significance is a more sound procedure than 
turning to a more restricted and less general model with 
the fewest number of fit parameters 22]. Specifically in 
our context, if the introduction of "f q as a fit parameter 
yields results expected from theoretical considerations, 
and at the same time this step raises statistical signifi- 
cance considerably, and furthermore makes the behavior 
of best fit parameters more consistent with expectation 
once experimental conditions are varied (collision energy 
and centrality than it is appropriate to conclude 

that current experimental data favor j q ^ 1. 

The methods of SH analysis are described elsewhere 
in great detail. We refer the reader to SHARE (statis- 
tical hadronization with resonances), the public SH suit 
of (FORTRAN) programs which we use in this analysis 
[25| . The parameter set of SHARE comprises also the fu- 
gacity Ai, which describes the asymmetry in the yield of 
up u and down d valance quarks. This parameters can- 
not be omitted (tacitly set to unity) in the AGS physics 
environment, since the initial state isospin asymmetry is 
not diluted by a very large produced particle yield. To 
understand the relevance of Ai, note that, in Boltzmann 
approximation prior to resonance decays (subscript 'o'), 
we have for the pion ratio: 

-?r = Aj = (K/X d ) 2 . (1) 

The reader will notice that we abbreviate particle yields, 
e.g., Y 7r ± by particle name tt . 

Using SHARE, we evaluate for a set of statistical pa- 
rameters the particle yields and as appropriate, ratios, 
and then find in least square minimization process the 
best parameter for the experimental results applicable to 
the considered collision system. We use here the most 
central event trigger available (5%). For the top energy 
at AGS, we adopt nearly the totality of the experimental 
results considered by prior authors as stated in Ref. [18J, 
listed here and comprising the following further develop- 
ments: 

i) Another way to introduce fugacity Ai is to conserve 
electric charge, see ^8| . Thus d > u asymmetry is equiv- 
alent to the requirement that the total charge Q of the 
hadronic fireball is a fixed fraction / = 0.391 of the to- 
tal baryon number B — B = 363 ± 10 (1) (we point out 
and count the number of measurements by having behind 
each value used in parenthesis a sequential number), es- 
tablished by proton and neutron content of the colliding 
gold ions, i.e., Q = 142 ±5 (2). We note that these values 
arise from our choice of the centrality trigger condition 
see, e.g., table I, in Ref. pr| . 

ii) We adopt a slightly different strategy in use the 



available particle yield data, in order to limit the in- 
fluence of the systematic error on the fit result: rather 
than to fit K~ = 3.76 ± 0.47 yield, we fit K+/KT = 
6.32 ± 0.65 (3) which is seen in Fig. 6 of Ref. p^ . 

iii) We explore what happens when we include in the 
fit the (redundant) but independently measured precise 
ratio K + /tt + = 0.197 ± 0.013 (4) (value taken from con- 
clusions to Ref. [27[ ) , which compares to the implied ratio 
K + /tt + — 0.177±0.037 derived from the individual yields 
we use: K+ = 23.7 ± 2.86 (5) (from table I in Ref. 0, 
including systematic error), and ir + = 133.7 ± 9.93 (6) 
(from extrapolation of Ref. [2^| as used in 01 ) The re- 
sults we present are obtained including this as an in- 
dependent measurement. Removing this ratio yields a 
reduction of P-value from, e.g., 65% to 55% but does not 
affect materially the results we discuss. 

iv) The other particle yield or, respectively, ratio values 
we use are: A = 18.1 ± 1.9 (7) , A = 0.017 ± 0.005 (8), 
and p/ir + = 1.23 ± 0.13 (9), in all of these following the 
procedure of Ref.0|. One of the reasons we did not 
deviate wherever possible from the data set of Ref . [Tsl] 
was to assure that the results we obtain can be compared 
directly with earlier work. 

v) We further study the influence of the recently pub- 
lished E917 ratio 0/K+ = 0.03 ± 0.006 [2<J, obtained 
within and interval of Sy = ±0.4 units of rapidity around 
the mid-rapidity. We choose from the other measure- 
ments presented (4> yield, (j)/ 7r ) this result as it is the 
most precise one, and suitable to a qualitative extrapo- 
lation to a total yield ratio. We assume that the Gaus- 
sian shape rapidity distribution has = 0.8 which is 
consistent with the data shown by E917 [2jj, and the en- 
ergy dependence systematics of obtained by the NA49 
collaboration, see Fig. 3 in [3(j ■ Since the interval of 
\5y\ = 0.4 < cr = 0.8 and <j = 0.8 < a K+ = 0.96 ± 0.6 
(cr K + is from Ref. |26j), to %-precision, the relation of 
yields within the rapidity interval to the Air total yield is 

Jl\ ~h±jL\ (2) 

K+ 1<* (Tk K+ l»s(-o.4,+o.4) 

Thus, we adopt as the 4tt ratio <j)/K + \ ilT = 0.025 ±0.006 
(10), which expresses the fact that the 4> is expected to 
have a narrower rapidity distribution than the K + . 

The SH model comprises aside of the already men- 
tioned 5 parameters T, /^b (equivalently valance light 
quark fugacity X q ), j s , j q , A/ also the fireball volume 
V and A s , the strange quark fugacity. In principle this 
last parameter can be fixed by the requirement that the 
strangeness and anti-strangeness content balances in the 
fireball, i.e., within the grand-canonical ensemble, 

(s) = (s). (3) 

However, we hesitate to use this constraint since: 

a) When chemical analysis includes both 7 S , j q and one 
of the sides of above equation involves in essence just one 
particle species (K + here) this constraint has in general 
several possible solutions. 
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FIG. 1: From top to bottom: statistical significance of the 
fit for the data set, see text, at a fixed value of 7 9 ; chem- 
ical freeze-out temperature T; baryochemical potential fiB', 
strangeness phase space occupancy divided by light quark oc- 
cupancy 7s/7<7, all as function of the prescribed to the fit 
light quark phase space occupancy f q . Dashed line in upper 
section includes in fit the <f>. Dotted lines are obtained impos- 
ing strangeness conservation. Vertical line, at j q = 1, is for 
orientation only. 



b) In the condition Eq. © , all the other measure- 
ments combine predicting the yield of K + required by 
the strangeness balance. This results in a consistency 
test of the data sets obtained from several experiments 
employing different instruments. 

All this speaks against implementing this constraint 
and the results we present were obtained not enforcing 
Eq. We will show by how much this constraint could 
be violated (in relative terms about 10%) and also show 
that its enforcement is not in essence altering the results 
shown. 

In this study, we have 7 free parameters and up to 10 
data points, of which one comprises partial redundancy. 
We will validate the fit instead of % 2 by the associated 
significance level P[%] which is obtained in SHARE us- 
ing the CERN library "PROB" procedure [p. P[%] 
is a function of the number of parameters p, and mea- 
surement points r. When these numbers are small, and 
in particular, the number of measurements is not much 



greater than the number of parameters, the value of to- 
tal x 2 must be much smaller than what is usually ex- 
pected for the 90% significance level result (\ 2 — p — r). 
Using this tool we find that 'good looking' figures arise 
for P ~ 15% applicable to the result shown in Ref. [l8[ . 
Such a low significance level could be result of chance, 
but more likely it means that either the error on some of 
the measurements considered is in realty bigger, or the 
theoretical model needs further refinement, such as, e.g., 
light quark chemical non-equilibrium. 

In view of the above, we first wish to understand if the 
7q = 1 choice (light quark chemical equilibrium) is com- 
pelling. Using the SHARE 1.2 package (without particle 
widths), we obtained the significance level (P- value) of 
our fits which we show in the top section of Fig.Q] We see 
(solid line) the significance level of our SHARE fit with 9 
first experimental results, which has a plateau near 80% 
in the range 0.3 < j q < 1.5. Eliminating the one redun- 
dant data point the significance level drops to 55% and 
the range of acceptable values of 7 9 widens further. The 
vertical line, in Fig. ^ is placed at the chemical equilib- 
rium, with assumed value j q = 1. 

The dashed line, in the top panel of Fig. ^ represents 
the significance level arising when the new 10th exper- 
imental point </>/K + is included in the analysis. The 
best fit is pushed well below chemical equilibrium with 
j q ~ 0.4, and the significance level rises nearly to 90%. 
As this result shows, there is no compelling reason to 
consider only the case of light quark chemical equilib- 
rium. Indeed this added experimental result favors chem- 
ical non-equilibrium for light quarks. The full chemical 
equilibrium j s — l,"f q — 1 appears inconsistent with the 
data set as can be seen inspecting the bottom panel of 
Fig. n where the ratio js/jq is presented. 

Dotted lines, in Fig.[I] were obtained enforcing 
strangeness conservation condition, Eq. ©, as was done 
in Ref. 18]. This condition does not influence in essence 
the discussion we present though in detail minor differ- 
ences arise. The two strategies we pursued in finding the 
best fits assure us that we have obtained the best fit in a 
situation involving many parameters, where several local 
minimal points are present. Our analysis shows that the 
suggestion that at low reaction energies the freeze-out of 
hadrons occurs well below the theoretical phase bound- 
ary between QGP and HG is not fully justified, as it is 
result of the assumption of chemical equilibrium. 

Specifically, we find for the best fit with <fi/K + a 
hadronization temperature T — 142 ± 3 MeV and the 
baryochemical potential /ib = 708 ± 60 MeV. We real- 
ize that both these values T and appear at first sight 
to be beyond the phase boundary region of the hadron 
gas, in the deconfined domain, given the estimate for the 
critical boundary [3^, and the related development of 
the liquid QGP phase model |3^. However, one has to 
realize that the relatively small values for j q and 7 S we 
obtained imply a significant reduction in the equilibrium 
QGP pressure, which pushes up the temperature of the 
non-equilibrium phase boundary. This opens the ques- 
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FIG. 2: Physical properties of the thermal fireball at chemical 
breakup as function of 7,. From top to bottom: thermal 
energy content weighted with -^/snn = 2.4 GeV; entropy per 
baryon S/ (B — B), strangeness per baryon s/(B — B) and the 
strangeness asymmetry of the fit, (s — s)/(s + s). 



tion if indeed the hadron yield results seen at AGS are 
not the product of the breakup of a deconfined state. In 
our earlier study of light ion collisions at AGS, we had 
noticed that the results allowed interpretation in terms 
of both confined and deconfined baryon rich fireball for- 
mation |34| . 

We can only speculate about the physical reasons be- 
hind these small values of favored phase space occupan- 
cies: 7 9 = 0.35 ± 0.27, 7 S = 0.23 ± 0.18. For example, if 
the deconfined state were to be well characterized at high 
baryon density by effective heavy quarks, with m ~ 350 
MeV, these quark equilibrium yields after hadronization 
at T = 142 MeV could well appear much below chemi- 
cal equilibrium. This is consistent with the low value of 
entropy per baryon, Fig. EI S/(B-B) = 14 ±2. This sit- 
uation is opposite to what happens for the small baryon 
density case (RHIC), with effectively small quark masses 
at hadronization, where the hadronization process over- 
populates the hadron yield. 

Considering together the AGS, SPS and RHIC energy 
ranges, we recognize that the hadronization temperature 
arises from a combination of several effects: the loca- 



tion of the equilibrium critical curve in the T, jUg-plane, 
shift in the critical curve due to chemical non-equilibrium 
conditions, and supercooling due to dynamics of the ex- 
panding fireball. The last effect is particularly signifi- 
cant at RHIC while the first is probably most relevant 
at AGS as we noted above. At RHIC, the fast trans- 
verse expansion is seen in many observables, and it is 
generally believed that a new partonic phase has been 
created. For an equilibrium system the phase cross-over 
would be expected at T = 164 ± 10 MeV How- 
ever, the wind of colored particles expanding against the 
color-non-conducting vacuum can go on until the sudden 
breakup near to T = 140 MeV |9|,|l5|- Aside of theoret- 
ical arguments, and fit results, such low hadronization 
temperatures are consistent with the observed resonance 
yield, for K\(f) H, as well as p°,/o and A(1520) 36]. 

Returning to the discussion of the physical proper- 
ties of the fireball: aside of entropy per baryon men- 
tioned above, we obtain other relevant physical prop- 
erties of the source of produced hadrons using the sta- 
tistical parameters to evaluate the phase space proper- 
ties. We present, in Fig. [21 from top to bottom, the 
thermal energy per baryons as fraction of the center of 
momentum energy available in the reaction, the single 
particle entropy per baryon, strangeness per baryon and 
the strangeness asymmetry in the fits we consider when 
strangeness conservation is not enforced. For the inter- 
esting range j q ~ 0.4, the value s/S ~ 0.008 ± 0.001 
is significantly below the SPS and RHIC level 0. This 
is consistent with the notion that strangeness yield rises 
with energy faster than light quark yield, as can be ar- 
gued considering the mass-energy threshold for the pro- 
duction of these flavors, and the time available in the 
collision. 

Considering that one measurement (0/K + ) pushes the 
discussion of AGS results toward a breakup of a decon- 
fined state, we wish to understand if other measurements 
could confirm and solidify this result. We present, in 
Figs. Eland 01 the variation of hadron ratios as function 
of 7 g with other relevant statistical parameters varying 
as indicated in Fig. We see considerable sensitivity of 
hadron ratios considered to the value of j q . We place the 
</>/K + experimental result in the top panel of Fig. El at 
an appropriate value of 7 9 , indicating with dotted lines 
the range of 7 9 consistent with the measurement error. 
Availability of other strange hadron ratios would provide 
the consistency check required to confirm that j q < 1. 

The results, presented in Figs. El and suggest which 
of the experimental ratios are of relevance in the study 
of chemical non-equilibrium. We see that some vary 
strongly, and that others vary little, e.g., A/p ratio is 
flat to within 10% in the entire range of j q considered. 
Clearly, experiments will have great difficulty to reach 
precision at %-level to see such variation. We considered 
also baryon to meson ratios involving S, A with K (and 
antiparticles) Mostly, these turned out to be very insen- 
sitive, i.e., flat to within 10%. The exception are ratios 
of A/K which are flat only for j q > 0.4. The variability 
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FIG. 3: From top to bottom: relative yields of 4>/K + , A/p, 
E~/A and fl/E~ as function of 7 9 . The experimental yield 
4>/K + is placed at best value of j q and the range (ls.d.) of 
possible 7, is indicated by vertical dots. 



for small 7 g is accounted for in the ratios shown in Fig.0] 
involving A. 

We note that the ratio of A/p ^1 at 7, < 0.2. 
The value of this ratio has been of considerable inter- 
est H3. The experimental result A/p = 3.6ttltll fa- 
vors a small value of j q We further note that the ra- 
tio 7r~/7r + , seen in Fig0J is consistent with the result 
7T-/7T+ = 1.23 ± 0.02 ± 15-20% obtained in central colli- 
sions > and again a small value of 7 g is favored when 
comparing the experimental result to the SH model pre- 
diction, obtained fitting other experimental data. 

This discussion confirms the importance of multi- 
strange hadrons and strange antibaryons as a valuable 
hadronic signature. In this context, it is important to 
understand the overall strangeness yield. Evaluating the 
number of valance strange quark pairs per baryon in 
terms of the parameters characterizing the phase space, 
we obtain for the entire relevant range: 



=- = 0.12 ± 0.02 ± 15%, 0.2 < 75 < 1. (4) 



(B-B) 

We show above the yield of s in emitted hadrons, since 
is greater than that of a in our study. When we enforce 



FIG. 4: From top to bottom: relative yields of 7r /ir + , A/p, 
H~/A and Q/S~ , as function of 7 q . 



strangeness conservation Eq. ©, the presented result is 
reproduced. The first error stated in Eq. (0} is the vari- 
ability of the yield as function of diverse parameters. We 
see that it is well below the estimate of the propagation of 
the experimental errors into the theoretical yield, which 
we present second as the estimated error. We note that 
ratio Eq. J3J includes aside of open strangeness, also a 
noticeable contribution contained in ss mesons {<j> and 
components in rj, 77'). 

We see, in Fig. |3 that there are more s-hadrons emit- 
ted than s-hadrons. If this s > s asymmetry were to be 
confirmed by more experimental data, this could be ex- 
plained by strangeness distillation phenomenon expected 
to occur at AGS energy j^. After the excess of s is 
evaporated, the residue is a quark soup enriched with s- 
quarks, a strangelett, which when metastable, may have 
escaped observation. 

For central collisions, at the top AGS energy, the trans- 
verse slopes of particle spectra T± ~ 200 MeV (e.g., K- 
spectra, see Ref. [2^ )■ This result leads to the average 
radial expansion of the AGS fireball at (v r ) — 0.5c pOj . 
Such a high speed of radial expansion is, for us, hard to 
understand. On the other hand, for a single freeze-out at 
T = 143 MeV we report here, we estimate the required 
(v T ) = 0.2-0. 3c (note that results of analysis do not 
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fully apply as these were done assuming chemical equilib- 
rium yields of resonances). The introduction of chemical 
non-equilibrium and the associated high chemical freeze- 
out temperature harbors the potential for consistent un- 
derstanding of particle spectra and yields, also within a 
kinetic reaction model [41j . 

In conclusion, we have presented a comprehensive and 
systematic study of the experimentally measured parti- 
cle yield ratios obtained in y / ijv7V = 4.8 Au-Au collisions 
(11.6 GeV/c Au on fixed target). We have analyzed the 
chemical freeze-out conditions allowing the QGP asso- 
ciated chemical non-equilibrium. We have found that 
the present data sample is not sufficient to argue deci- 
sively for or against the QGP presence at AGS energy 
scale. However, consideration of the recently measured 
yields of <p, the tt~ /tt + along with indirectly evaluated 
A/p favors as the result of the data analysis the chemi- 



cal non-equilibrium hadronization at the phase boundary 
between the confined and deconfined baryon rich phase. 
The results we presented indicate that exploration of the 
phase transition between baryon rich confined and de- 
confined phases may be possible by means of relativistic 
heavy ion beams in the energy range of 10 A GeV/c. 
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